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ABSTRACT 



Timmes. Brown. & TruranI found that metallicity variations could theoretically ac- 
count for a 25% variation in the mass of ^^Ni synthesized in Type la supernovae (SNe 
la), and thus account for a large fraction of the scatter in observed SN la luminosities. 
Higher-metallicity progenitors are more neutron-rich, producing more stable burning 
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products relative to radioactive ^^Ni. We develop a new method for estimating bolo- 
metric luminosity and ^^Ni yield in SNe la and use it to test the theory with data from 
the Supernova Legacy Survey. We find that the average ^^Ni yield does drop in SNe 
la from high metallicity environments, but the theory can only account for 7%-10% 
of the dispersion in SN la ^^Ni mass, and thus luminosity. This is because the effect is 
dominant at metallicities significantly above solar, whereas we find that SN hosts have 
predominantly subsolar or only moderately above-solar metallicities. We also show 
that allowing for changes in 0/Fe wit h the metallicity [Fe/H1 doe s not have a major 



effect on the theoretical prediction of iTimmes. Brown. & Truranl so long as one is 



using the 0/H as the independent variable. Age may have a greater effect than metal- 
licity - we find that the luminosity weighted age of the host galaxy is correlated with 
^^Ni yield, and thus more massive progenitors give rise to more luminous explosions. 
This is hard to understand if most SNe la explod e when the prirn aries reach the Chan- 



drasekhar mass. Finally, we test the findings of [Gallagher et al.L that the residuals of 
SNe la from the Hubble diagram are correlated with host galaxy metallicity, and we 
find no such correlation. 

Subject headings: supernovae: general — surveys — galaxies: abundances — galax- 
ies: highredshift — galaxies: fundamental parameters 



1. Introduction 



Type la supernov ae (SNe la) make excellent cosmological standard candles (iRiess et al.lll998 : 



Perlmutter et a 



shape (|Phillips 



19991), though o nly after they are empirically ca librated based on their lightcurve 



19931) and color (|Riess. Press. & Kirshneilll996|) . There is consensus that the lu 



minosity of T ype la supernovae (SNe la) arises from the radioactive dec ay of ^^Ni to ^^Co, and 
finally to ^^Fe ( Truran. Amett. & Cameron|[l967 : Colgate & McKeeiri969 ). Therefore, the primary 
variable controlling SN luminosity is the amount o f ^^Ni synthesized during the explosion, though 
as second-order effects t he location of the ^^Ni (e.g. Branch & Khokhlovll 19951: Nugent etaP 19971; 
Pinto & Eastman llOO Oa^. the o pacity of the overlying material (e.g. Kh okhlov. Mueller. & Hoeflich 



1993 



Pinto & Eastman 



the ejecta (|Howell et al 



2000bl: 



20011 : 



Mazzali et al. 



Kasen et al 



200ll : lKasen & Woosley. 2007) . and asymmetries in 



20041) can redi stribute the energy in wavelength, time . 



or space. However, because SNe la are largely spherical dWang et al.lll996l : [Leonard et al.ll2005l) . 
the bulk of ^^Ni is located near the c enter, and simple a ssumptions about the opacity are a rea- 
sonable approximation, Amett's Rule (|Amettil 197911 19821) can be used to e mpirically estimate the 
^ ^Ni mass frorn the o ptical-infrared lightcurve to an accuracy of 5 - 15% dStritzinger et al.ll2006l : 
Blinnikov et allbooeh . 



-3- 



The variables that control the yield of ^^Ni in a SN la are not well understood, although 
several factors are known to be important. Iron-peak elements are produced when burning occurs 
in a state of nuclear statistical equilibrium, which happens at the highest densities near the center 
of the white dwarf. If this burning occurs in a neutron-rich region, more neutron-rich, stable, iron- 
peak elements like ^^Ni and ^^Fe will be produced relative to ^^Ni. Since main-sequence stars with 
higher m etallicity abundances produce whit e dwarfs with more of the neutron-rich nuclides ^^Ne 
and ^^Fe, iTimmes. Brown. & TniranI (120031 hereafter TBT03) hypothesized that high metallicity 
progenitors will produce less ^^Ni and therefore less luminous SNe la. 

In the absence of any identified SN la progenitor, the best one could hope for is to see traces 
of the progenitor's metallicity in the unburned outer laye rs of the SN by looking at Type la spec- 
tra from only a f ew days after the explosion in the UV faoflich. Wheeler. & Thielemann| [l998: 
Lentz et al.ll200lb . Unfortunately the difficulty of obtaining UV spectra has prevented the gath- 
ering of a statistically significant sample of spectra from the earli est phases. In the case of local 
SNe la it has not been possible to schedule HST quickly enough (|Foley. Filippenko. & Jhall2008 : 
Sauer et al.ll2008b . Another option is to observe higher redshift SNe la in the restframe UV, though 
their faintness at such early epochs rel ative to their host galaxy makes observations earlier than one 
week before maximum light difficult (lEUis et al.ll2008|) . More importantl y, even when UV spectra 
have been obtained their theoretical interpretation is not straightforward (|Ellis et al.ll2008h . 



Since a true progenitor metallicity has so far been unobtainable, the next best thing is to 
determine the average metallicity in the environment of the SN. The two most commonly used 



appro aches hav e been to measure spe ctroscopic line indices (|Hamuy et al.ll2000l : [Gallagher et al 



20081) or ratios (iGallagher et al.ll2005h to determine a global average host galaxy metallicity, or to 
take advantage of the fact that galaxies o ften have metallicity gradients, and study projected galac- 
tocentric offset as a proxy for metallicity dwang. Hoflich. & Wheeleill 19971 : llvanov. Hamuv. & Pintd 
2000; IGallagher et al.ll2005h . The drawbacks of the spectroscopic line techniques are that they re- 
quire complicated corrections for contamination by emission, and that the interpretation of the line 
ratios as met allicity requires mapping observations to a sometimes-incongruoii s grid of theoretical 



models (e.g. iTrager et al.ll2000l : 



Kewlev & Dopitall2002l : IGallagher et all 120081) . Furthermore, the 



line ratios only represent the current average metallicity - the metallicity when the SN la progen- 
itor formed may have been different. In addition, high signal-to-noise ratio spectra are required, 
and important lines like Ha are redshifted into the IR at z > 0.5. These factors limit the pos- 
sible sample size of suc h a data set. The drawbacks of the offset technique a re that projection 



effects can be confusing dWang. Hoflich. & Wheeledl 19971 : 



van den Berghlll997h. as is differ ential 



reddening between SN e near and far from galaxy center s (|Hatano. Branch. & Deatonlll998l) . and 
Malmquist-like effects (jHowell. Wang. & Wheeled 120001) . Further complications include that dif- 
ferent galaxy types have different metallicity gradients, and the current location of the SN may not 
be representative of where the progenitor formed. 
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With the exception of lGallagher et al.1 (|2008|) . hereafter G08, none of the above studies found a 
significant correlation between metallicity and SN la properties, although it is not clear whether this 
is caused by the lack of an effect or whether the techniques used so far do not have the required level 
of precision. G08 compared absorption-line strengths of 29 E/SO galaxies which hosted SNe la to 
a grid of models and found that galaxies with high iron abundances host less luminous supemovae. 
They went on to conclude that age played a greater role, since galaxies with characteristic ages 
greater than 5 Gyr produce ~ 1 mag fainter SNe la in V. Their most controversial finding is that 
the MLCS lightcurve fitter could not fully correct for SN luminosities in different metallicity envi- 
ronments; when Hubble residuals were correlated with metallicity, they found a significant trend. 
But since these findings are based on a small sample size, rely on an uncertain metallicity determi- 
nation techniques, and rely entirely on a lightcurve fitter (MLCS2k2: |jha. Riess. & Kir shner 2007) 
that has recently shown spurious systematic trends with Hubble residuals ( iconley et al.ll2007|) . they 
are worth scrutinizing with independent techniques and a much larger sample size. 

Here we use an independent approach not previously applied to testing SN luminosity-metallicity 
relations — we estimate the host galaxy metallici ty from its mass using spectral energy distribution 
(SED) fits to the host galaxy photometry and the lTremonti et al.l (|2004l) mass-metallicity relation. 



Another discriminator of the metallicity effect is that it may cause lur ninosity differences a s 
a function of redshift, since cosmic metallicity is increasing with time (e.g. ISavaglio et al.ll2005|) . 
This is in contrast to many other possible contributors to SN la luminosity (for example, white 
dwarf core density) which would not be expected to change with redshift. We perform the first 
tests of the effect of inferred metallicity on SN la nickel mass with high redshift supernovae, using 
data from the Supernova Legacy Survey (SNLS). 

In this paper we first update the TBT03 SN la metallicity-^^Ni mass relation to account for 
the fact that 0/Fe can vary as a function of Fe/H (§2). Then we test the theory using data from 
the SNLS. To do this, in §3.1 we estimate host galaxy masses from PEGASE-2 model fits to 
UMgMrMiMZM (subscripts denote Megacam filters) host galaxy photometry. In §3.2 we use galaxy 
mass-metallicity relations to estimate host metallicities. In §3.3 we determine SN ^^Ni mass from 
SNLS photometry. Finally, we test the theory with the data in §4 and discuss the implications in 
§5. 



2. The Variation of ^°Ni Mass with Progenitor Metallicity 

TBT03 note that a white dwarf's ^^Ne abundance should trace its progenitor star's O abun- 
dance since C, N, and O are converted to ^^N during hydrogen burning (^'^N(j9, 7) is the slowest 
step in the CNO cycle), and this ^^N is converted to ^^Ne during core helium burning. Since ^^Ne 
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is neutron-rich, its presence during the burning in the Type la supernova favors the production of 
stable, neutron-rich nuclides like ^^Ni and ^"^Fe at the expense of radioactive ^^Ni, the decay of 
which powers the supernova. Thus, more metal-rich progenitors are expected to produce dimmer 
SNe la. 

The progenitor star's abundance of ^^Fe is also important, but one possible limitation of 
TBT03 is that metallicity was treated as a single variable, i.e. TBT03 make the simplifying as- 
sumption that 0/Fe is constant relative to Fe/H. In reality 0/F e varies as a function of Fe/H (wit h 
a great deal of scatter), at least in the solar neighborhood (e. g. lWheeler. Sneden. & Truranlll989|) . 
and the degree of variation is expected to vary with the history of star formation. Here we refine 
the estimate of the mass of ^^Ni to include changes in the relative abundances of O and Fe. 

TBT03 showed that M55, the mass of ^^Ni synthesized in the explosion, depended on the 
electron abundance Ye, 

M56 = Ml(,{5SYe-2S). (1) 

Here M^^ is the mass of ^^Ni synthesized at an electron abundance Yg = 1/2, i. e. for a pure C-0 
white dwarf. It does not include material that has been neutronized by in situ electron captures at 
high densities. To a good approximation, we can consider the white dwarf to consist of ^^C, ^^O, 
^^Ne, and ^^Fe. The electron abundance is then 



Ye = 1(1- 22yNe - 56yFe) + lOFNe + 26Yf. 



6) 



(2) 



where 7, = Xi/Aj is the abundance of nuclide i, Xj and A, being the mass fraction and atomic 
mass number, respectively. Since ^^Ne is derived from the aboriginal CNO abundances of the 
progenitor main-sequence star, = + + ^^c- Nitrogen is a sec ondary element, and we 
therefore write Fn = >b [1 + (N/0)], with (N/0) = 0.0316 + 126(0/H) dvila Costas & Edmunds 



19931). This formula includes the secondary pr oduction of nitroge n at high metallicity and agrees 
with derived abundance ratios from the SDSS (iLiang et al.ll2006|) . Carbon is a primary element, 
but is mostly produced in lower-mass stars, which tends to m ake the ratio (C/Fe) nearly con- 
stant at subsolar metallicities, a t least down to [Fe/H] > —2 (|Chiappini. Romano. & Matteucci 
2OO3I : iBensby & FeltzingI 120061) , and there is evidence th at [C/Fe] decreases slightly at super- 
solar metallicity in the thin di sk fee nsby & Feltzin^2006 ). We take (C/Fe) = 8.7, its solar value 
([Asplund. Grevesse. & Sauvalii2005i) . 



Writing Yq =Zh(0/H) and Fpe =^H(Fe/H), settingZn = 0.7392 dAsplund. Grevesse. & Sauval 



20051), and inserting equation ^ into equation ([T]) gives 



1 - 58 < 0.763 



+ 93.1 



+ 7.92 




(3) 
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The factor of 58 comes from the assumption that the two dominant nuclides in the bum to nuclear 
statistical equilibrium are ^^Ni and ^^Ni. In actuality ^"^Fe is also present, and the ratio of ^'^Fe 
to ^^Ni depends on the thermal conditions of the explosion (TBT03), so that the coefficient varies 
between 58 and 54. As shown by TBT03, this results is a small scatter (< 1% at Z = 3Zq) about 
the relation in eq. ([3]). 

The two abundance ratios in equation ([3]) are correlated, albeit with a great deal of scat- 
ter. We can therefore eliminate Fe/H in favor of 0/H and the coefficients from the relation be- 
tween [Fe/H] (the "meta llicity") and [0/Fe]. Here we are us i ng the standard notation [A/B] = 



log(A/B) — log(A/B)0. iRamfrez. AUende Prieto. & Lambert! (|2007[) fit a set of linear relations, 
[0/Fe] — a -l-Z? [Fe/H], for the thin disk, the thick disk, and the halo (see Table [T]). We therefore 
write Fe/H = (Fe/H)Q lOl^'^/^l, use the linear relation between [0/Fe] and [Fe/H] and the identity 
[O /Fe] = [0 /H] - [Fe/H] to obtain [Fe/H] ^ ( [O /H] - a) / ( 1 + Z?) , and insert the solar photospheric 
abundances dAsplund. Grevesse. & Sauv3l2005l) . log (O/H)0 = -3.34 and log (Fe/H)^ = -4.55, 
into equation ([3]) to obtain 



M56 



1-0.044 



(0/H) 



10 



-3 



1+0.122 



(0/H) 

10-3 



+ 



(0.19+0.53/7+fl)/(l+i) 


r(o/H)i 






10-3 





(4) 

For many purposes, it is more convenient to work directly with [Fe/H] as the primary variable. We 
therefore write (0/H) = {O /H) qIO^'^/^^^^^^^/^^ and use the linear relation between [0/Fe] and 
[Fe/H] to obtain 

^ = 1 -0.020 X lO^+^i+^^P'^/"] |l +0.056 X io«+(i+*)P«/h] +0.64 x io---h[Fe/n]\ (5-^ 

In both equations, (Hj) and ([5]), setting a = b = will give the case of an syste m where O /Fe is held 

at its s olar value. Figure[T]shows the equation (HI) with coefficients taken froml Ramfrez. AUende Prieto. & Lambert 
(l2007h for the thin disk {dashed line) and thick disk {dash-dotted line). We also show, for compar- 
ison, the relation when 0/Fe is fixed at its s olar value {solid line). This last curv e differs slightly 
from that in TBT03 because we are using the lAsplund. Grevesse. & Sauvall (l2005h abundances for 
the solar composition, and we include the secondary production of ^"^N. This scaling is appropriate 
for the Milky Way, but supernova hosts may have a different metallicity history. This does not 
appear to be much of a concern, however, since the differences between the predictions using the 
thick and thin disk parameters are small for [0/H] < 0.5. We conclude from this section that 0/Fe 
variations can safely be neglected at the level of precision of the TBT03 theory, for the relevant 
range, [0/H] < 0.5. We note from Figure[i|that TBT03 would have predicted a variation of nickel 
mass with metallicity of less than 10% had they restricted their considerations to the observation- 
ally relevant range of [0/H] < 0.5 that corresponds to [Fe/H] < 0.6 in the thin disk model. 
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Fig. 1. — Theoretically expected variation in M56, the mass of ^^Ni synthesized, as a function of 
0/H. The solid line shows the prediction of TBT03, which fi xed [0/Fe] — 0. The other lines show 
the pr ediction when 0/Fe varies according to linear relations iRamfrez. AUende Prieto. & Lambert 
(|2007i) fitted to the thin (dashed lines) and thick (dot-dashed lines) disk. For illustrative purposes, 
we show each over a two decade span in metallicity as plotted in TBT03, though the thick and thin 
disk relations were derived over narrower metallicity ranges (see Table [U, and SN hosts do not 
have average metallicities spanning this full range (section |4]). 
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3. Method 



3.1. Estimating host galaxy masses 

To test the TBT03 theory, we must determine host galaxy metallicities, and to do that we 



outlined in 


Sullivan et al. 


(2006 


), we used the code Z-PEG ( 


Le Borgne & Rocca-Volmerangell2002) 


to fit PEGASE.2 ( 


Fioc & Rocca-Volmerange 


1997, 


1999 


) models to host galaxy UMgurMiMZM 



photometry from the SN 
analysis of ISullivan et al. 



( Sullivan et al 



2008 



.S to determine their masses and star formation rates. We extend the 



(l2006h to include data from th e SNLS 3rd year cosomology sample 



Conlev et al.ll2008bl : lOuv et alJl2008h . 



The SNLS revisits 4 one-degree fields five times per lunation for the 5-6 lunations per year that 
each field is visible. Images with seeing < 0.8" and no SN light are combined to provide a deep, 
multi-year reference image from which the host galaxy properties are measured. SN types and 



host galaxy redshifts were determined from Gemini, VLT, and K eck spectroscopy (I Howell et al 



20051 : iBronder et al.ll2008l : iBaumont et al.ll2008l : EUis et al.ll2008h . We categorize each galaxy by 
its star formation rate per unit mass (specific star formation rate; sSFR) into strong star formers 
(log(sSFR) > —9.5), weak star formers (—12 > log(sSFR) > —9.5), and passive galaxies with no 
measurable star formation rate (log(sSFR) < —12). Here, SFR is the average star formation rate 
over the last 0.5 Gyr, which is estimated from the optical SED fitting. We assign errors correspond- 
ing to the range of masses that can be fit within a model given the photor netric errors, or the sc atter 
between acceptable model fits in terms of whichever is larger. As in lSuUivan et al.l (|2006h . we 
only use SNe at z < 0.75 where the SNLS sample is relatively complete. For more inforrnation o n 
the SNLS, including discoveries, the photometric system, and SN typing, see lAstier et al.l (|2006|) . 



3.2. Estimating metallicity 



To determine average host galaxy metallicities we use the iTremonti et al.l (|2004h correlation 
between host galaxy stellar mass and gas phase 0/H as determined from 53,000 galaxies in the 
Sloan Digital Sky Survey. There is a 0. 1 dex scat t er in th is relationship. For galaxies with log < 
8.4 we use the low mass extension of 
for zero redshift 



Lee et al. J2OO6 I). Both of these relationships were derived 



Savaglio et al.l (|2005|) find that higher-redshift galaxies h ave lower avera ge gas 



phase metallicities. They give an "effective Mass" versus r e dshift relat ion (iSavaglio et al.l (12005) 
equation 9), which we use here to adapt the Tremonti et al. (2004) and Lee et al. (|2006 ) relations 
to the redshifts of our galaxies. 



Technically the 



Tremonti et al, 



(|2004|) relationship only covers gas phase metallicity, and is 
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not guaranteed to be applicable in elliptical galaxies, where there is generally little cold gas. Still, 
there are several reasons not to dismiss the sample of passive galaxies from this analysis. First, 
the amount of star formation in these galaxies is inferred from their broadband colors, not from 
spectroscopy or morphology, so the presence of gas is not ruled out. Second, because the mech- 
anism for retaining metals is the depth of the galaxy gravitational potential well, one might still 
expect ellipticals to follow this trend, independent of the presence of gas. The SN progenitor was 
bom when there was star formation, and thus much more gas in the galaxy. Finally, we caution 
that stellar metallicity and gas phase metallicity are not identical, though they are correlated (e.g. 



Cid Femandes et al.ll2005 ). 



3.3. Estimating ^^Ni mass 



The mass of 5^Ni can b e estiniated from th e peak bolometric luminosity and the rise time of 
a SN la using Amett's Rule (|Amettill979l . ll982h . that the luminosity radiated at maximum light is 
approximately equal to the energy deposited by radioactive decay: 

^bol 



Mm 



rs{tR) ' 



where y is the ratio of bolometric to radioactiv ity luminosities, near unity (e.g. Branch 



1992; 



Hoflich & Khokhlovl 19961: Blinnikov et al l2006h . Here we adopt a value of 1 . 2^=0.1 JNue ent et aP 
1995, : .Branch & Khokhlovll 19951^ to be consistent with previous SNLS work (IHowell et al.,.20061 



though a global change to this number does not affect our conclusions. S is the radioactivity lumi- 
nosity per solar mass of ^^Ni from its decay to ^^Co and subsequent decay to ^^Fe: 



7.74 X lO^^e"'^/^-' 



+ 1.43 X 1043[g-Wiii _g-W8.8] M 



© 



Here tR is the time in days for the supernova to rise from explosion to maximum light 
tR — sx 19.5 days (IRiess et al.lll999l : IConley et al.l 120061) . where s is the " stretch", a scaling fac 
tor for the time axis of the SN lightcurve that correlates with luminosity (jPerlmutter et al.lll997 



IConley et al. 2008a ). Note that simply taking the mass difference between ^^Ni and ^^Co, as some 
authors have done, will overestimate the energy production. The heating fro m the decay o f ^^Ni 
is primarily from the radiative decay of the 1.7 MeV nuclear level in ^^Co (|jundelll999[) . We 
include contributions from the pos itron decay branching in the heating from the decay of ^^Co 
Jcolgate. Petschek. & Kries3ll980l) . 



The least straightforward part of determining ^^Ni mass from SN la lightcurves is estimating 
the bolometric luminosity from incomplete data. For rare, very well-observed supernovae with data 
spanning much of the UVOIR (ultraviolet, optical, infrared) region, all of the flux can be added, and 
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only a small correctio n is necessary for missing data (e.g. IContardo. Leibundgut. & Vaccall2000l : 
Stritzinger et al.ll2006|) . However, this has to date only been possible for small numbers (~ 15) of 
low-redshift supemovae. Because the best-observed supemovae tend to be the closest and bright- 
est, and thus not in the Hubble flow, these SNe often have uncertain distances and luminosities. 

Other authors have used the lumino sity in a single band, combined wi th a bolometric correc- 
tion, to estimate the bolometric flux (e.g. lBranchlll992l : iHowell et al.ll2006|) . This allows the study 
of greater numbers of supemovae, achieving enhanced global accuracy, at the expense of individ- 
ual accuracy. It also makes possible the study of high redshift supemovae where the restframe NIR 
is inaccessible. 

Bolometric corrections work well at maximum light for Type la supemovae, because the peak 
of emission is always in the restframe B band, near 4000A (in the absence of strong reddening). 
Figure 121 (bottom p anel) shows a typical SN la spectral energy distribution at maximum light (from 
Hsiao et al. ( 2007 1). updated with NIR data). The upper panel shows the cumulative flux as a 
function of wavelength. Only ~ 5% of the SN la flux is emitted in the UV shorward of the U filter, 
and less than 10% is emitted in the IR0 Since roughly one-quarter of the SN la flux is emitted in 
the B band, and this fraction is relatively constant from SN to SN, it is easy to see why the B band 
is the most commonly used indicator of SN la luminosity. 

Here we adopt a new approach to bolometric flux estimation that exploits the excellent multi- 
wavelength coverage of SNLS supemovae. For each SNLS SN la, we smoothly adjust (i.e. "warp") 
the Hsiao SED template so that it has the observed fluxes in gMfMiMZM- We then integrate the 
template SED over the wavelength interval where we have observations, and correct for the unob- 
served fraction. We turn this bolometric observed flux into a luminosity using the redshift and a 
flat cosmology with Hq=10 km s^^ Mpc^^ and = 0.3. An example is shown in Figure [31 



Note that this process uses much of the same machinery as the SiFTO lightcurve fitter (|Conley et al 



2008a() . The first step, warping the SED to match the observed colors, is identical. But in the 
next step the methods diverge - SiFTO calculates the restframe flux of SN through another filter, 
whereas here we calculate the integrated flux over the region spanning all observed filters. 

Mathematically, we define / as the fraction of the bolometric flux covered by our observations, 
such that: 

^ jH{X)dX ^ ^ 

where H{X) is the Hsiao SN la SED template, X\ is the wavelength at 1/5 the filter height of the 



'Our template spans lOOOA — 2.4 /xm - since there is no evidence for significant SN la flux outside of this range, 
we assume it to be neghgible. 



0.001 



5000 10000 15000 20000 25000 

Wavelength [A] 



Fig. 2. — Bottom: The lHsiao et al.l (|2007h SN la maximum-light spectral energy distribution tem- 
plate (plotted on a log axis), updated to include NIR data. Top: Cumulative flux as a function of 
wavelength. Johnson-Cousins filters are shown with labels corresponding to the percentage of total 
flux through each. Approximately 3/4 of the flux for a typical SN la is radiated in the optical. 
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Restframe Wavelength [A] 
2000 3000 4000 5000 6000 7000 




2000 4000 6000 8000 ^ 10000 12000 
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Fig. 3. — Determination o f the bolometric flu x for an example SN la, SNLS-03Dldt at z = 0.612. 
The dashed line shows the lHsiao et al.l(|2007|) SED, and the solid line shows it warped to match the 
observed gMfMiMZM fluxes. This warped SED is then integrated over the observed region (shaded), 
from Ai to A2, arbitrarily defined to start and stop where the bluest and reddest filters are 1/5 of 
their maximum height. Since this region contains only 74% of the total Hsiao SED bolometric flux 
(see Fig. 12]), the integrated flux is divided by 0.74 to obtain the bolometric flux. 
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bluest filter for which there is data, and A2 is the wavelength at 1/5 the filter height of the reddest 
filter for which there is data. The choice of 1/5 of the filter height as the cutoff is arbitrary - 
changing it does not affect our results. 

We define 5(A) as the SED after it is warped to match the observed colors, so that the bolo- 
metric flux is 



Fboi 



f 



(7) 



Flux errors are given by the SiFTO covariance matrix (see lConley et al. Errors on the other 

quantities are: 



CJ56 



Ni 



1.4 X 1044^-2.2. _ 2.5 X 1042^-0.18., ^ 



752 



(8) 



(9) 



Based on constraints from theory (e.g. lBranchlll992l: iHoflich & Khokhlovlll996l) . and from esti- 
mating ^^Ni using different empirical methods ( Stritzinger et al.l 2006 ). we adopt ay = 0.1 and 
generally this term dominates the error budget. 

In the absence of dust extinction, one should use the raw absolute magnitude (i.e. not corrected 
for stretch or color), because this reflects the true luminosity of the SN. The complication is that the 
observed color is a mixture of reddening due to dust (which should be corrected for), and intrinsic 
SN color (which should not be corrected). If the color correctio n were primarily due to normal dus t, 
it would be expecte d to follow the M i lky Way extinction law ( Cardelli. Clavton. & MathisI 19891) 



but it does not (e.g. ^ Wang et al.ll2008l : lElias-Rosa et al.ll2008l : iKrisciunas et al 



200 



Conlev et al. 



2007!) . There are the oretical reasons to expect, and sorne evidence, that the rnajority of SNe la suf- 
fer little extinct ion (iH atano. Branch. & Deatonl 1 19981 : IPerlmutter et al.l 1 19991 : IPhillips et aP [l999: 
Comminsll2004l) . though some of the reddest o utliers surely have si gnificant dust extinction. As our 
primary result, we chose to make a color cut dSuUivan et al.ll2008l) to eliminat e the most extincted 
SNe rather than making an uncertain extinction correction Jconley et al. 2007 ). The parameters of 
the color cut are given in Tabled The chosen cut parameters, and whether or not the cuts are made 
at all, do not affect the conclusions. We also study the effect of making a color correction, which 
corrects for dust, but also improperly corrects intrinsic luminosity dispersion. Such results give a 
lower limit on the dispersion of ^^Ni. 

To make the color correction we followed the sa me steps as above, but first we use the 



color excess measured for each SN la to unredden the Hsiao et al 



(120071) S ED template using 



the lCardelli. Clavton. & MathisI (119891) law with Ry = \.% (iConlev et al.ll2007D. We use a fiducia l 
color of c = 0, where c is a linear combination of restframe B — V a nd U —B Jconley et al. 2008a ). 
roughly corresponding io E{B — V)q — —0.057 (lAstier et al.ll2006|) . 
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Table 1. [0/Fe] for different populations 





a 


b 


range 


Thin disk 


0.096 ±0.004 


-0.327 ±0.016 


-1.0 < [Fe/H] < -0.3 


Thick disk 


0.370 ±0.027 


-0.121 ±0.043 


-0.8 < [Fe/H] < 0.3 


Halo 


0.388 ±0.049 


-0.048 ±0.071 


-1.4 < [Fe/H] < -0.4 


Note. 


— Coefficients 


of the fit fO/Fel 


= a + blFe/H] from 



Ramirez. AUende Prieto. & Lambert) (|2007|) 



Table 2. Parameters of color cut 





color 


a 


b 


c 


U02-B 


0.1856 


-0.0003 


-0.1681 


0.0790 


U-B 


0.4031 


0.1816 


-0.0567 


0.0477 


V-R 


1.0129 


0.1957 


-0.5634 


0.0646 



Note. — SNe la rejected if the B — V color is 
redder than a (color +0.5) +b{s— 1) +c, where 
s is the stretch (lightcurve width). "?702" is a 



virtua l filter useful in k-correcting (I Hsiao et al. 



2007|) high reds hift data, defined to be bluer 



than the U filter dSuUivan et al.ll2008l) . AB-V 
and at least one other color was required for all 
SNe la. Cuts were made in as many filters as 
data were available. The cuts mimic those used 



in ( Sullivan etal.ll2008 ) 
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Mathematically, if Hq{X) is the unreddened template, then this process can be summed up by 
the following two equations: 

f = ^-h L_ (10) 

and 

Ftoi= y ■ (11) 



4. Results 



4.1. Galaxy Mass vs. SN la stretch 



It is well know n that SN lightcurve widths c orrelate wit h host galaxy morphology (IHamuv et al 

19951 12OOO : HoweiboOll : loallagher et aDl2005l). ho st color jBranch. Romanishin. & Baronll 19961) . 
or specific star formation rate ( Sullivan et al.ll2006 ). in the sense that SNe with wider lightcurves 
are more commonly associated with late-type galaxies and vice versa. The fact that elliptical 
galaxies are generally larger (in a spatial sense) than spirals also may account for trends seen 



between lightcurve w idth and projected galactocentric distance (|Wang. Hoflich. & Wheeleiill997 



Gallagher et a l."2005l). And while SN rate s have been studied as a function of host galaxy mass 
([Mannucci et al.ll2005uSullivan et al.ll2006|) . no study has examined lightcurve width as a function 
of host galaxy mass. Figure |4] shows such a comparison for SNLS SNe la. Keeping in mind that 
passive galaxies tend to have a high mass, the results are as expected from previous work — high 
mass galaxies host few high stretch SNe la, and low mass galaxies host few low-stretch SNe la. 



4.2. ^^Ni Mass vs host galaxy metallicity 

In Figure [51 we have converted galaxy mass to gas phase metallicity and converted integrated 
flux to 5^Ni mass as outlined in ^|3l We have overplotted the TBT03 relationship between metallic- 
ity and predicted SN la ^^Ni yield, corrected for 0/Fe differences as outlined in ^ In this figure 
we make a color cut, but not a color correction. Figure |5] shows that there is a drop in average ^^Ni 
mass for SNe la from high mass, high metallicity (12-i-log(0/H)> 8.8) galaxies. TBT03 predicted 
a 25% difference in ^^Ni yield over a factor of 3 difference in progenitor metallicity. However, 
most of the effect occurred at metallicities 12-i-log(0/H) > 9.2 (see Fig.[T]), whereas most SNe la 
in this study occur in galaxies with lower metallicity. Figure [5] shows that the TBT03 prediction 
is fairly flat over the range of metallicities determined for actual SN hosts. Even when the steeper 
"thin disk" relation is used from Fig. [H at most, the expected effect is only 0.06 Mq in this range. 
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Log (Host Mass) 

Fig. 4. — SN relative lightcurve width (stretch) as a function of host galaxy mass determined from 
fitting PEGASE.2 galaxy models to host galaxy UMgurM^MZM photometry. Circles (red online) 
have hosts with no detectable star formation (log sSFR < — 12), where sSFR is specific star forma- 
tion rate. Star symbols (blue) denote strong star formation (log sSFR> —9.5), and squares (green) 
are in between. 
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Fig. 5.— 5^Ni mass as a function of inferred average host gas phase metallicity. The dashed line is 
the prediction of TBT03 modified to use the thin disk variation in 0/Fe as described in ^J2l Points 
are SNLS SNe la with ^^Ni mass and inferred average host metallicity determined as described in 
the text and colored as in Fig. HI Vertical error bars include flux errors, stretch errors, and errors 
associated with Amett's Rule. Horizontal error bars include errors fro m the model fits, but for 
clarity do not include the 0.1 dex scatter from the lTremonti et al.l (120041) relation ship. The gap at 
12+log(0/H)=8.3 ari ses from the disc ontinuity between the lTremonti et al. (2004) mass-metallicity 
relationship and the Lee et al.l (120061) low mass extension (there is no gap in host galaxy masses 
as can be seen from Fig. |4]). Points are plotted by metallicity, and are corrected for redshift mass- 
metallicity effects as mentioned in the text, but the upper axis gives an approximate corresponding 
mass at zero redshift. Triangles are averages in 0.2 dex bins. Outer error bars are the standard 
deviation of the points in that bin, and the inner error bars are the errors on the mean. 
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In comparison SNe la show a wide range in derived ^^ Ni mass, here 0.2 to 1.0 M(7^, though ranging 
from 0.1 to 1 3Mn for extreme cases like SN 199 Ibg Jpilippenko et alJI 19921) and SNLS-03D3bb 
jHowell et alJbooel) . 



In other words, the TBT03 theory appears to be qualitatively consistent with observations, 
though there is additional scatter not explained by the theory. An additional error of 0.16 M© of 



56t. 



Ni had to be added to the Ni error bars in quadrature to achieve J /DOF = 1. 

Several SNe la with ^^Ni below 0.2 Mq were conservatively eliminated with the color cut, 
on the grounds that extinction might be causing a low derived ^^Ni mas s, although is also poss ible 
that these are subluminous SNe, which have intrinsically redder colors (|Gamavich et al.ll2004h . 



Various effects can systematically alter the scaling of the observed ^^Ni mass in Fig. |5l includ- 
ing changing the Hubble constant Hq or the ratio of radiated to deposited radioactive decay energy 
7, but these do not affect the overall results, because the normalization of the TBT03 prediction is 
arbitrary. Here we have chosen the normalization by fitting to the data. 

Another way to visualize the difference in ^^Ni yield between high and low metallicity galax- 
ies is the histogram in Figure[6l Supernovae in more metal-rich hosts (12-i-log(0/H)> 8.8) have a 
much lower average ^^Ni yield than those in metal-poor hosts. The high metallicity set has a peak 
in the distribution near 0.4 Mq, while the low metallicity set peaks near 0.7 Mq. The probability 
that the ^^Ni yields from high and low metallicity hosts were drawn from the same distribution, 
according to the Kolmogorov-Smimov (KS) test is 9 x lO^'^. 



4.3. Reddening correction 

It is possible that the above results overestimate the scatter in SN la ^^Ni mass, if significant 
dust extinction remains after the color cut described in ^ To test this, we make the color correction 
in Figure |71 As expected, the scatter is reduced - now only 0.12 Mq must be added in quadrature 
to the error bars to give ;f^/DOF = 1. Since this color correction includes an intrinsic SN la color 
correction as well as a dust correction, it will over-correct the SN la magnitud es, but we do it to ge t 
an estimate of the limiting case. If SNe la are color-corrected using /3 = 2.8 ('C onlev et allboOTl) . 
the scatter is reduced to a range of 0.6 Mq in ^^Ni for all but a few outliers. Thus the upper limit 
on the TBT03 metallicity effect is 0.06/0.6 = 10%. 
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Fig. 6. — Histogram of SN •'"Ni mass split by inferred host gas phase metallicity. Higher metallic- 
ity hosts have a lower ^^Ni yield, on average. This uses the non-color-corrected data from Fig. \5\ 
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Fig. 7. — Same as Fig. [H but a color correction has been made. This correction would be appro- 
priate if all SNe la were the same color, and all reddening and dimness were due to dust. However, 
there is an intrinsic color-luminosity relationship for SNe la that it is not currently possible to sep- 
arate distinctly. Thus a color correction overcorrects the luminosity, and the result is a lower limit 
on the scatter in SN la ^^Ni mass. 
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Fig. 8.— ^^Ni mass versus host luminosity weighted age as estimated from PEGASE.2 models. A 
color cut on the SN magnitudes has been made, but not a color correction. The line shows a fit to 
the data with slope —0.15 ±0.03 and correlation coefficient -0.38. 
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Fig. 9. — Same as Fig. [8l but a color correction has been made. The line shows a fit to the data 
with slope —0. 1 1 ± 0.028 and correlation coefficient -0.37. 
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4.4. Progenitor Age 



Since there is an age-metallicity degeneracy, it is difficult to separate effects that may be 
caused by metallicity from those that may be caused by age of the progenitor. An age-brightness 
effect must exist for SNe la, since both SN rate and average luminosity are incr eased in young stel- 
lar populations (e.g. lHowellll200ll : iMannucci et al.ll2005l : ISuUivan et al.ll2006[) . Figure [8] confirms 
that there is a correlation between SN ^^Ni mass and host luminosity-weighted age as estimated 
from PEGASE.2 models — the slope of the line is significant at ~ 5a. This correlation is un- 
derstandably noisy - the PEGASE.2 luminosity-weighted age is merely an indicator of the mean 
age of an entire galaxy - the SN progenitor may have a very different age. Any age-brightness 
relation must be related to with the mass of the secondary star - in an old population more massive 
secondaries do not exist. But this is difficult to understand in the Chandrasekhar mass model - if 
all SNe la explode when the primary nears the Chandrasekhar mass, why should the mass of the 
secondary matter? This might be an indication that the binary evolution and mass transfer history, 
which is strongly dependent on the mass of the secondary, influences the outcome of the explosion. 



4.5. ^^Ni mass vs. metallicity grouped by stretch 



Figure [To] shows ^^Ni mass vs. 12-i-log(0/H) in three stretch bins. The most striking finding is 
that low stretch (s < 0.9) SNe la almost exclusively occur in galaxies with high implied metallicity, 
roughly above-solar. At the highest stretches (s > 1.05), the inferred ^^Ni mass is relatively insen- 
sitive to inferred host galaxy metallicity. Also, the scatter in ^^Ni mass is lower for low stretch SNe 
la. This is likely due to several factors. For the low-stretch SNe, on the high side the ^^Ni mass is 
limited by the stretch cut, and on the low side it is limited by selection effects - it is difficult to find 
and spectroscopically confirm the lowest-stretch SNe. But there is also some evidence that SNe 
la in passive galaxie s, which predominantly host low stretch SNe la, make better standard candles 
(|Sullivan et al.ll2003r) . Some of this lower dispersion is likely due to lower dust extinction in these 
hosts, but it is also possible that mod erately low-stretch SNe la are intrinsically more uniform as a 
population than high stretch SNe la (|Howell & Nugentll2004l) . 



4.6. Hubble residuals 

G08 found a trend between SN la host galaxy metallicity and residuals from the Hubble di- 
agram, indicating the lightcurve fitter used, MLCS2k2, cannot correct for the full ra nge of SN la 



prope rties. Since the demographics of supemovae are known to evolve with redshift ( Howell et al 



20071), any problem fully correcting SN la magnitudes could have dire consequences for SN la cos- 
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Fig. 10. — Ni mass versus implied 0/H as in Figure [51 split by stretch. Within a stretch bin there 
is no obvious inferred host galaxy metallicity dependency (line shows TBT03 prediction). The 
bottom panel shows that low stretch SNe la are always produced in higher-metallicity galaxies, 
and the scatter in ^^Ni mass is lower for low-stretch SNe la. 
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Fig. 1 1 ■ - — Hubble residua l vs. inferred host galaxy metallicity for SNe la from the SNLS 1st year 
sample (lAstier et alj|2006l) . The solid line is a fit to the data, though it is consistent with zero slope 
at the ~ la level. The dotted line shows t he G08 relation , here ruled out at > 99.9% confidence. 
The metallicity dispersion of 0.1 dex in the lTremonti et al.l (|2004r) relation was added in quadrature 
when calculating the fit, but is not plotted for clarity. Residual errors do not include intrinsic SN la 
scatter about the Hubble line. No redshift or color cut was applied to the data; applying either or 
both gives a slope consistent with zero and still rules out the G08 relation at > 99.9% confidence. 
Symbol shapes and colors are the same as Fig. HI The upper axis shows the log mass for a galaxy 
with the metallicity found on the lower axis at z = 0. The convention used is the same as that in 
G08 — overluminous SNe have a negative Hubble residual. Note that we are using an indicator of 
[0/H] as the independent variable, whereas Gallagher use an indicator of [Fe/H]. Any systematic 
trend with Hubble residual should show a slope with respect to either variable. 
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mology. [Gallagher et al.1 cite a 9% systematic error in the meas urement of the Dark Energy equa 



tionof state, w, larger than the current best statistical errors on w (|Kowalski et al. lAstier et al 
20061) . 



However, using our independent metallicity indicator and the SiFTO lightcurve fitter we find 
no systematic trend in Hubble residual vs. inferred host galaxy metallic ity. The results are shown 
in Figure [TTl where the SNe la include those from lAstier et al.l (l2006b sample, but we apply no 
redshift or color cut. The solid li ne shows the best fit to the data, based on a Markov Chain Monte 
Carlo analysis (LINMIX in IDL; iKellyl 120071) . which takes into account errors on the x-axis, and 
is allowed to add extra scatter to the data (here 0.1 1 mag) to achieve the best fit. The slope of the 
best-fit line, —0.10 ±0.07 is consistent with at the ~ la level, and is inconsistent with the G08 
slope of 0.26 at 5a. In 10,000 Monte Carlo simulations, none produced a slope as steep at that 
found by G08. 

Possible explanations for the difference between our results and those of G08 include that 
SiFTO can correct for differences in SN la properties that the current version of MLCS2k2 cannot, 
that one or the other metallicity indicator produces systematically flawed results, or that there are 
selection effects such that we are not comparing similar samples. We examine each hypothesis in 
turn. 



SiFTO vs. MLCS: Recently Ijha. Riess. & Kirshneil (120071) found evidence for a "Hubble 
Bubble" — local SNe la fit with MLCS2k2 give a different Hubble constant than do SNe la well 
into the Hub ble flow. IWood-Vasev et al.l (|2007[) considered this to be a 6% systematic error on 
w. However. IConley et al.l (l2007b did not see the same effect with either the SiFTO or SALT2 
(|Guy et al.ll2007b lightcurve fitters, and traced the discrepancy to problems in the way MLCS2k2 
handled dust extinction. MLCS2k2 assumed color excess beyond any assumed intrinsic SN la color 
relation is due to i^v = 3.1 (average Milky Way dust), though this appears not to be appropriate 
(SALT2 and SiFTO do an empirical fit to determine color-luminosity relationships). GOB do not 
provide enough information to determine whether a spurious reddening correction is driving their 
result. 

Flawed metallicity indicators: Since we cannot make a definitive statement about whether 
the j-axis of Figure \TT\ (or the equivalent figure in G08) is flawed, the next culprit to examine 
is the X-axis. No metallicity indicator is perfect — the various caveats throughout the text attest 
to the limitations of our adopted method. Estimating a metallicity through a mass-metallicity 
relationship and determining the mass from photometry is certainly indirect. Still, it should be 
correct on average, especially since we are probing several orders of magnitude in mass, ranging 
from< 10^ to> 10i\and more than 1 dex in metallicity. 



A stronger argument is that there seems to be no way to rearrange the points on the x-axis of 
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Figure[TT]to produce the G08 relation. On the other hand, it is possible to imagine that the x-axis of 
the upper left panel of G08 Fig. 9 could be systematically flawed. The four points with the highest 
metallicity drive the relation. These four SNe all have E/SO hosts with ages less than 3 Gyr, among 
the lowest in the G08 sample. Thus, the hosts of these SNe la are the most likely to have H/3 
emission contaminating the H/3 absorption line index. Nominally this is used to determine age, but 
because of the age-metallicity degeneracy, it also affects the metallicity determination. As can be 
seen from GOB Figure 2, when the emission correction is applied 1/3 to 1/2 of their sample either 
falls off or is moved onto their theoretical age-metallicity grid. Even after correction, eight of their 
SNe la are more than la too high off the grid, and their metallicities are extrapolations. 

Selection effects: Since G08 target nearby SNe la in E/SO hosts, they study many SNe la 
too dim to show up in a magnitude-limited sample like the SNLS. But the four SNe with high 
metallicities driving the G08 Hubble residual trend are all of normal magnitude — 19 < My < 
— 19.5, not the type absent from the SNLS sample. The metallicity range is also similar for the 
two studies, spanning a range ± ~ 0.5 dex relative to solar. But if the differences between the two 
studies were due to some unknown selection effect, an important point is that it will not produce 
a systematic effect on cosmology, because the SNLS sample reflects cosmology as it is actually 
practiced, and we have shown that it produces unbiased Hubble residuals. 



5. Conclusions 



We find that host galaxies with a higher inferred metallicity produce SNe la with less ^^Ni 
on average. The effect proposed by TBT03, increased neutronization in higher metallicity envi- 
ronments, can qualitatively explain this trend, though there remains additional scatter in SN la 
luminosity that cannot be explained by the theory. TBT03 predict that SNe la in high metallic- 
ity environments should be 25% dimmer than those in low metallicity galaxies, though most of 
the dynamic range in this prediction occurs at metallicities several times solar. We have shown 
that SNLS supernovae occur in galaxies where the implied average metallicity is a few tenths of 
a dex above solar or below, and in this range the theory predicts only a 0.06 M0 change in ^^Ni 
mass. The average ^^Ni yield does appear to decrease in metal-rich galaxies in accordance with 
the theory, but the large scatter (0.2 - 1.0 M©, or ~ 0.6 Mq with color correction), indicates that 
the theorized metallicity effect has at most a 7%-10% effect on ^^Ni yield. 

The luminosity of S Ne la appears to be dominated by other effects, possibly those caused by 
the ag e of the progenitor dHowelibooTI : iGallagher etaDl2005l : iMannucci etaPbood : ISuUivan et aL 



20061) . We find a correlation between host galaxy luminosity-weighted age and SN ^^Ni yield. 
This indicates that more massive stars (explicitly, massive secondaries, since they set the timescale 
for explosion in binaries), produce more luminous explosions, a result not well understood if most 
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SNe la explode near the Chandrasekhar mass. We have also shown that allowing for changes 
in 0/Fe with the metallicity [Fe/H] does not have a major effect on the theoretical prediction of 
TBT03 for metallicities within the observed range of our sample, so long as one is using the 0/H 
as the independent variable. 

Unlike G08, we find no trends between residuals from the Hubble diagram and host metal- 
licity. We conclude that metallicity is not a significant systematic error for the measurement of 
w. The GOB finding may be due to problems with their lightcurve fitter, MLCS2k2, or with their 
method of host metallicity determination which relies on an uncertain correction for H/3 emission, 
and the extrapolation of theoretical models to cover the host galaxies of SNe la which fall off the 
grid. 

The se findings are b r oadly consistent with the range of local SN la host galaxy metallicities 
found by lOallagher et al.l (|2005|) . though we differ in how we compare the model to the data, and 
in the interpretation of the results. Despite using completely different methodology and data, both 
studies find that the majority of SN hosts have a metallicity between 8.4 < log (0/H) + 12 < 
9.1. Additionally we find several lower metallicity, low mass hosts. These are unde rrepresented 
in low redshift searches which target large, higher-mass galaxies (e .g. iLi et al.ll200lb . But rather 



than convert the observations to ^^Ni mass, [Gallagher et al. 



prediction into an expected difference in Amis (|Phillipslll993|) . They found a fairly steep evolution 



( 2005 ) converted the TBT03 model 



in the Amis model prediction over the observed range in host 0/H, in contrast to the flat prediction 
over the same range found here (Fig. [5] ) and in TBT03. This is beca use they used th e dela yed 
detonation models of Hoflich et al. (2002) and the empi rical relations of Gamavich et al.l (|200 4) to 
convert the model to the observe d parameter Ami s . The Hoflich et al. ( 2002 1) models had a varying 
detonation transition density, so [Gallagher et al.l (|2005|) found a steeper relation than one would 
derive if only 0/H were varied while holding the other properties of the explosion fixed. As a 
result, they concluded that the TBT03 prediction is a poor fit to the data, whereas we find that it 
could account for some of the variation in SN la properties, but it is not the dominant effect. 



Our findings are also consistent with lPiro & BildstenI (l2008f) . who roughly estimated that the 
TBT03 theory could on ly explain a small fra ction of SN la luminosity based on average SNLS SN 
la stretches rep orted in iHowell et al.l (|2007r) and average SN host metallicity estimated from the 
galaxy types in lSuUivan et al. hood) . 



One caveat regarding our findings is that we are using a proxy for the global average galaxy 
metallicity - we have no way of knowing the actual progenitor's metallicity. One could imagine 
an effect that would skew the results, such as an enhanced SN rate in metal-poor or metal-rich 
environments, or similar diminution of the rate in either environment, such that the average SN 
metallicity would not be representative of the average metallicity in the g alaxy. While there has 



been speculation from time to time that such effects could in principle exist (IKobayashi et al.lll998 
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Langer et alj|2000|) . there is no observational evidence that they do. 



The absolute calibration of the mass-metallicity relationship is also uncertain (|Kewley & Ellison 



20081), but changing it would have the main effect of moving all points in Fig. [5] to slightly lower 
or higher metallicities. Our conclusions would remain unchanged - that the metallicity effect can 
only account for a fraction of the scatter in SN la luminosities. 



Piro & BildstenI (|2008l) find that weak reactions during low-level burning in the white dwarf 



in the ~ 1000 years leading up to explosion (i.e. simmering) may increase the neutron excess. 
This effect may set a "floor" to the level of n eutronization that dorn inates over the TBT03 metal- 



licity effect at low metallicities Z/Z© < 2/3. IChamulak et al.l (|2008|) also find that simmering can 



increase neutronization. Since neutronization effects during simmering depend on the degree of 
simmering, and at least to first order do not appear to depend on metallicity, it is possible that sim- 
mering could wash out the TBT03 metallicity effect. Significant uncertainties, including the role 
of the convective Urea process, remain, and further study is required. At the very least it seems 
that simmering cannot enhan ce the TBT03 metallic ity effect enough to account for the significant 
scatter in SN la luminosities (IPiro & Bildstenll2008h . 



One additional caveat is that it is possible that a change in metallicity could alter the explosion 
mechanism itself, not just whether the process produces more ^^Ni or ^^Ni. IChamulak. Brown. & Timmes 
(|2007|) found that laminar flame speed also depends on the metallicity, but this is probably a higher- 
order effect. At any rate, this process is so poorly understood that it is not yet possible to test this 
effect with observations. 
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